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Measurement of the Inclusive Jet Cross Section in pp Interactions at = 1.96 TeV 
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We present a measurement of the inclusive jet cross section in pp interactions at y/s = 1.96 TeV 
using 385 pb _1 of data collected with the CDF II detector at the Fermilab Tevatron. The results are 
obtained using an improved cone-based jet algorithm (Midpoint). The data cover the jet transverse 
momentum range from 61 to 620 GeV/c, extending the reach by almost 150 GeV/c compared with 
previous measurements at the Tevatron. The results are in good agreement with next-to-leading 
order perturbative QCD predictions using the CTEQ6.1M parton distribution functions. 

PACS numbers: 13.87.Ce, 12.38.Qk, 13.85.Ni 


The differential jet production cross section at the 
Tevatron probes the world’s highest momentum trans¬ 
fers in particle collisions, is potentially sensitive to a 
wide variety of new physics, such as quark composite¬ 
ness Q, and tests perturbative QCD (pQCD) over more 
than eight orders of magnitude. There was great inter¬ 


est when the inclusive jet cross section measured by the 
CDF collaboration at the center of mass energy y/s = 1.8 
TeV 111 exhibited an excess in the high transverse en¬ 
ergy Ep region when compared to next-to-leading order 
(NLO) QCD predictions obtained using then-current par- 
ton distribution functions (PDFs) Q. Global PDF anal- 


4 


ysis by CTEQ 3 E] demonstrated that the excess could 
be explained by an enhanced gluon distribution at high 
momentum fraction x (x > 0.3). Recent global PDF fits 
(CTEQ6, CTEQ6.1, MRST2004) 011 @ , which include 
the Run I Tevatron jet data '2J, UJ|> find an increased 
gluon density at high x and provide a good description 
of the Run I Tevatron data. The gluon distribution is 
still poorly constrained at high x (see e.g. Ref. (3) and 
contributes significantly to the theoretical uncertainty for 
many interesting processes at the Tevatron and the LHC. 
The increase in \Js from 1.8 to 1.96 TeV, together with 
higher luminosity in Run II, allows more precise jet pro¬ 
duction measurements with a significantly extended kine¬ 
matic range. 

Jet algorithms cluster together objects such as par- 
tons or particles, or energies measured in calorimeter 
cells. The clustering algorithms rely on the association 
of these objects based on transverse momenta (the kx al¬ 
gorithm), or angles (the cone algorithm), relative to a jet 
axis. A measurement using the kx algorithm is reported 
in Ref. 0 - In this letter, we report the results of an 
inclusive jet measurement usin g a cone algorithm for the 
rapidity region 0.1 < |t/| < 0.7 |l2j. Cone jet algorithms, 
rather than kx algorithms, have been used dominantly at 
hadron collider experiments mainly due to the simplicity 
in constructing corrections for the underlying event and 
for multiple interactions in the same bunch crossing 0 
It is worth noting that, previously, results from a cone 
algorithm 0 and kx algorithm :14j by the D0 collabo¬ 
ration showed only marginal agreement at low px where 
corrections for multiple interactions and underlying event 
become important. We use the Midpoint algorithm, an 
improved iterative cone clustering algorithm 0 . It is 
difficult to use previous iterative cone algorithms 110 
with higher order pQCD calculations due to the presence 
of infrared singularities. The Midpoint algorithm places 
additional seeds between stable cones having a separa¬ 
tion of less than twice the size of the clustering cones; 
the use of these additional seeds reduces the problem 
with infrared singularities. 

The CDF II detector is a magnetic spectrometer which 
is described in detail elsewhere 0 . Here we describe 
briefly those components that are crucial to this measure¬ 
ment. The central detector consists of a silicon vertex 
detector inside a cylindrical drift chamber. Surround¬ 
ing the tracking detectors is a superconducting solenoid 
which provides a 1.4 T magnetic field. Outside the 
solenoid is the central calorimeter, covering a pseudo¬ 
rapidity (77) 0 range up to 1.1. The central calorimeter 
consists of 48 modules, segmented into towers of granu¬ 
larity Ar) x A cf> « 0.1 x 0.26 and divided into electromag¬ 
netic (CEM) and hadronic (CHA) sections. The CEM 
is a lead-scintillator calorimeter; the CHA is an iron- 
scintillator calorimeter with a depth of approximately 4.7 
interaction lengths. The energy resolution of the CEM 
for electrons is a(Ex)/Ex = 13.5%/i/£r(GeV) © 2% 


while the average energy resolution of the CHA for 
charged pions is a(Ex)/Ex = 50%/ \JEx(GeV) © 3%. 
The forward region, 1.1 < |7y| < 3.6, is covered by 
the “Plug Calorimeters” consisting of lead-scintillator for 
the electromagnetic section and iron-scintillator for the 
hadronic section. The region between the central and 
forward calorimeters, 0.7 < |7y| < 1.3, is covered by an 
iron-scintillator hadron calorimeter with similar segmen¬ 
tation to the central calorimeter. 

This measurement uses a data sample corresponding to 
an integrated luminosity of 385 pb _1 collected between 
February 2002 and August 2004. The data were collected 
using four trigger paths. The Level 1 trigger requires a 
calorimeter trigger tower, consisting of two calorimeter 
towers adjacent in 77, to have either Ex > 5 GeV or Ex > 
10 GeV. At Level 2, the calorimeter towers are clustered 
using a nearest neighbor algorithm. Four trigger paths 
with cluster Ex > 15, 40, 60, and 90 GeV are used. 
Events in these paths are required to pass jet Ex > 20 
(J20), 50 (J50), 70 (J70), and 100 (J100) GeV thresholds 
at Level 3, where the clustering is performed using a cone 
algorithm with a cone radius R cone = 0.7. 

Cosmic ray events are rejected by a cut on the miss¬ 
ing transverse energy (Jfix) significance 0|. For J20, J50, 
J70, and J100, we remove events having a/?T significance 
greater than 4, 5, 5, and 6 GeV 1 / 2 , respectively. In addi¬ 
tion, all events containing jets withpr > 360 GeV/c and 
passing the analysis cuts were visually scanned, and no 
cosmic ray events were found. The efficiency of the fix 
significance cut is 100% for jets at low px (65 GeV/c) and 
decreases to 92 % for jets at high px (550 GeV/c). We 
reconstruct ^-vertices by fits to tracks in the event and a 
beamline constraint, and select the vertex with the high¬ 
est total px of the associated tracks as the event vertex. 
In order to ensure that the jet energy is well measured, 
the event vertex is required to be within 60 cm of the 
center of the detector along the beamline. The efficiency 
for this cut is determined to be 95% from the beam pro¬ 
file measured using a minimum bias sample. Jets are 
required to have a rapidity |j/| between 0.1 and 0.7 to re¬ 
duce the effects of the gap between calorimeter modules 
and at the transition region between the central and plug 
calorimeters. 

There are two essential stages for any jet algorithm. 
First, the objects (partons, particles, or calorimeter tow¬ 
ers) belonging to a cluster are identified. With the Mid¬ 
point algorithm the cluster is a cone of radius 0.7 in (77, <f>) 
space. For reasons dealing with problems of unclustered 
energy endemic to iterative cone algorithms 0 , the clus¬ 
tering radius is at first set to R con e/ 2(= 0.35), and then 
later expanded to its full size as discussed below. Sec¬ 
ond, the kinematic variables defining the jet are calcu¬ 
lated from the objects comprising a cluster. The Mid¬ 
point algorithm makes use of four-vectors throughout 
the clustering process. The four-vector for each tower 
is computed as a sum of vectors for the electromagnetic 
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and hadronic compartments of the tower; the vector for 
each compartment is defined by assigning a mass-less vec¬ 
tor with magnitude equal to the deposited energy and 
with direction from the event vertex to the center of each 
compartment 01- The detector towers are sorted in or¬ 
der of descending pr- Only towers passing a seed cut, 
ptouier > pf^ ed ^ are used as starting points for the initial 
jet cones. The seed threshold is chosen to be 1 GeV/c; 
its value has a negligible effect on jets in the kinematic 
region used in this measurement. A tower passing the 
threshold of 100 MeV/c is clustered into a cone and even¬ 
tually into a jet if the separation from the axis of the cone 
in (p, cp) is smaller than R cone /2. There is no threshold 
for particle and parton clustering. After each iteration 
the jet centroid position is updated. The jet clustering 
is repeated until all of the jet cones are stable. A cone is 
stable when the tower list is unchanged from the previous 
iteration. After all stable cones have been determined, 
the clustering radius is expanded to the full size ( Rcone )• 
The use of the smaller initial cone results in an expected 
cross section approximately 5% larger due to the inclu¬ 
sion of jet energy that would have remained unclustered 
in the default Midpoint algorithm (0. At this point, an 
additional seed is defined at the midpoint between any 
two cones separated by less than 2 R cone and the iteration 
process is repeated. Two overlapping cones, if present, 
are merged into a single jet if the shared energy is larger 
than 75% ( fmerge = 0.75) of the energy of the jet with 
lower pt; otherwise the shared towers are assigned to the 
nearest jet. This splitting/merging procedure is iterated 
until the tower assignments to jets are stable. The jet 
kinematic properties are defined using a four-vector re¬ 
combination scheme Q. 

The inclusive differential jet cross section is defined as: 

d 2 a _ 1 1 N jet /e 

dprdy Ay f Ldt A pt 

where Nj et is the number of jets in the pt range Apr, 
e is the trigger, $t significance cut and z -vertex cut ef¬ 
ficiency, f Ldt is the effective integrated luminosity, and 
Ay = 1.2 is the rapidity interval used in the analysis. 
A trigger efficiency greater than 99.5% is required to in¬ 
clude the jets collected by a given trigger threshold. The 
measured calorimeter level jet cross section must be cor¬ 
rected for detector effects and for energy from additional 
pp interactions in the same bunch crossing (pile-up). For 
this sample, the average number of additional pp inter¬ 
actions is about 0.9. The pile-up corrections subtract 
0.93(±0.14) GeV/c for each additional 2 -vertex from the 
measured jet pt [l8| . 

The detector response corrections are determined from 
a detector simulation and a jet fragmentation model. The 
detector response is determined using a GEANT-based de¬ 
tector simulation 19] in which a parametrized shower 
simulation (gflash |20|) is used for the calorimeter sim¬ 
ulation. The GFLASH parameters are tuned to test-beam 


data for electrons and high-p^ charged pions and to the 
collision data for low-p^ charged hadrons 0 . PYTHIA 
6.216 0, with Tune A jEj). I23j. is used for the pro¬ 
duction and fragmentation of the jets. Tune A refers to 
the values of the parameters describing multiple-parton 
interactions and initial state radiation which have been 
adjusted to reproduce the energy observed in the region 
transverse to the leading jet in the data from Run I. It 
has also been shown to provide a reasonable description 
of the measured energy distribution inside a jet |2i|. 

The measured pr spectrum must be corrected for de¬ 
tector effects before it can be compared to theoretical 
predictions. We cluster the final state stable particles 0] 
in PYTHIA using the same algorithm as the one used to 
cluster calorimeter towers. The resulting jets contain 
all the particles from the pp collision, including those 
from the hard scatter, multiple parton-parton interac¬ 
tions and beam remnants. The correction, done in two 
correlated steps, is determined from a large sample of 
PYTHIA events, passed through the CDF detector simu¬ 
lation. First, a pT-dependent correction is determined by 
matching the particle jet to the corresponding calorime¬ 
ter jet and is applied to each measured jet. A binned 
spectrum is formed from the corrected pt of each jet. The 
bin widths are chosen commensurate with jet energy res¬ 
olution and statistics. The pt correction ranges from 1.17 
at low p T (65 GeV/c) to 1.04 at high p T (550 GeV/c). 
The spectrum must be further corrected for bin-to-bin 
jet migration due to the finite energy resolution of the 
calorimeter. This unfolding correction depends on the 
detector energy resolution and the true spectrum as well 
as the pT-dependent correction that was applied in the 
first step. The PYTHIA events are reweighted to match 
the experimental spectrum before the correction factors 
are calculated. A bin-by-bin unfolding correction is then 
determined by taking the ratio of the binned hadron level 
cross section and calorimeter level cross section corrected 
by the pr-dependent correction. The size of the unfolding 
correction varies from 1.30 at low pt to 2.31 at high pr- 
The applied corrections remove the detector effects from 
the raw cross section and the corrected hadron level cross 
section can now be compared to theoretical predictions. 

The main systematic uncertainties on the measured 
inclusive jet cross section arise from four sources: the 
jet energy scale, the jet energy resolution, the unfolding 
of the measured cross section to the hadron level cross 
section, and the luminosity. The dominant source of un¬ 
certainty is from the jet energy scale. The energy scale is 
known to better than 3% over the entire transverse mo¬ 
mentum range, leading to an uncertainty on the jet cross 
section varying from 10% at low pt to —39% at high p T , 
comparable to the uncertainty achieved by CDF in Run 
I. The uncertainty due to the jet pr resolution is deter¬ 
mined by the pt resolution difference between the data 
and the PYTHIA Monte Carlo. The uncertainty on the 
cross section varies from about 6% at low pt to about 
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10% at high pt . The uncertainty associated with the un¬ 
folding correction is determined by correcting a herwig 
6.5 12511 dijet sample using the corrections derived from 
the PYTHIA sample. This uncertainty is determined to 
be less than 5% at high px and lesss than 10% at low 
Pt- The luminosity uncertainty is 6%, independent of 
Pt and is not included in the quoted systematic error. 
Other effects considered were determined to have a neg¬ 
ligible effect on the cross section. Adding all of these 
contributions in quadrature yields a total experimental 
systematic uncertainty on the inclusive jet cross section 
varying from approximately 15% at low px to approxi¬ 
mately % at high pt- 

To compare the data with predictions for jets of par- 
tons as obtained from NLO calculations, the data must 
be further corrected for underlying event and hadroniza- 
tion effects. It is also possible to correct the NLO predic¬ 
tions for the same effects; the two corrections are simply 
the inverse of each other. For the former, we correct for 
the energy in the jet cone not associated with the hard 
scatter, i.e., from the collisions of other partons in the 
proton and antiproton. The latter corrects for particles 
outside the jet cone originating from partons whose tra¬ 
jectories lie inside the jet cone. It does not correct for the 
effects of hard gluon emission outside the jet cone, which 
are already accounted for in the NLO prediction. The 
bin-by-bin hadron-to-parton corrections are obtained by 
applying the Midpoint clustering algorithm to the hadron 
level and to the parton level outputs of the PYTHIA Tune 
A dijet Monte Carlo samples, generated with and without 
an underlying event. The sample without the underlying 
event was generated by turning off multiple parton in¬ 
teractions. The underlying event correction results in a 
decrease of the cross section varying from 22% at low px 
to 4% at high px\ the hadronization correction increases 
the cross section by 13% at low pr, and by 3.5% at high 
Pt- herwig provides consistent results on the hadroniza¬ 
tion corrections, but predicts smaller underlying event 
energy; the difference in the underlying event correction 
is taken as the underlying event correction uncertainty. 
In previous measurements at the Tevatron UM , the 
hadronization corrections were not applied to the data. 
The inclusive jet cross section is shown in Fig. ^ and 
Table □ lists the cross sections with the statistical and 
systematic uncertainties at the hadron and parton levels. 
Also included in Tableware the explicit factors applied to 
the hadron level cross section to obtain the parton level 
cross section. The experimental and theoretical jet cross 
sections are obtained by averaging over the transverse 
momentum bins. 

Current NLO theoretical predictions for inclusive jet 
production exist only at the parton level, for which the 
final state consists only of 2 or 3 partons [H, m m 
For our comparisons with theory we use the calculation 
of EKS (2^|. The ratio of the inclusive jet cross section, 
corrected to the parton level, to the NLO QCD predic¬ 


tions using the CTEQ6.1M PDFs is shown in Fig. [2] The 
Midpoint jet algorithm has been applied to the 2 or 3 par- 
tons in the final state of the EKS calculation. In order 
to mimic the properties of the splitting/merging step, 
present at the experimental level but not at the NLO 
parton level, a parameter R sep with a value of 1.3, has 
been introduced E2 Two partons are clustered within 
the same jet if (1) they are within R con e (0.7 for this mea¬ 
surement) of the jet centroid and (2) within R sep x R CO ne 
of each other. The use of R sep = 1.3 results in a reduc¬ 
tion of the theoretical cross section prediction by approx¬ 
imately 5%, roughly independent of jet transverse mo¬ 
mentum, as compared to the prediction obtained when 
Rsep is not used in the calculation. In the EKS predic¬ 
tions, the renormalization and factorization scales (pr 
and pf) have both been set to p!^ 1 / 2. Using a scale of 
p^ (2 p^ 1 ) rather than p?* /2 leads to a theoretical pre¬ 
diction for the jet cross section lower by approximately 
10% (20%) over the entire pt range and a larger x 2 in 
the global PDF fits [t| . The gluon distribution has been 
determined in the global fits, primarily by the Tevatron 
Run I jet data, using a renormalization and factoriza¬ 
tion scale of Px*/ 2. Thus, for self-consistency, this scale 
should be used in the NLO comparisons. 

We show in Fig. [5] the experimental uncertainties 
for the inclusive jet cross section and the theoretical 
uncertainties estimated from the 40 CTEQ6.1M error 
PDFs [ 2 ] • The PDF uncertainty is the dominant the¬ 
oretical uncertainty for most of the transverse momen¬ 
tum range. The correction for underlying event and 
hadronization is model dependent. The error associated 
with this correction is added in quadrature to the to¬ 
tal experimental error and shown in Fig. |2]as the outer 
shaded band. The data are in good agreement with the 
NLO QCD predictions, which is consistent with what is 
reported in Ref. & 

It is important to emphasize that the CTEQ6.1M 
gluon density is already “enhanced” at high x and so au¬ 
tomatically leads to a larger prediction for the jet cross 
section than older PDFs such as CTEQ5M. Also shown 
in Fig. Q is the prediction using the latest PDF set from 
the MRST group The MRST2004 PDFs also con¬ 
tain an enhanced higher x gluon, leading to reasonable 
agreement with the CDF jet measurement. 

In conclusion, we have measured the inclusive jet cross 
section in the range 61 < px < 620 GeV/c using an 
improved iterative cone clustering algorithm, Midpoint. 
The new measurement extends the jet transverse momen¬ 
tum range over previous measurements at the Tevatron 
by about 150 GeV/c. The data are well described by 
NLO QCD predictions using CTEQ6.1M PDFs, within 
the theoretical (PDF) and experimental uncertainties. 
No new physics is indicated in the high px region. Inclu¬ 
sion of these data in future global PDF fits will provide 
further constraints on the gluon distribution at large x. 

We thank the Fermilab staff and the technical staffs 


7 


of the participating institutions for their vital contribu¬ 
tions. This work was supported by the U.S. Department 
of Energy and National Science Foundation; the Italian 
Istituto Nazionale di Fisica Nucleare; the Ministry of 
Education, Culture, Sports, Science and Technology of 
Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the 
Republic of China; the Swiss National Science Founda¬ 
tion; the A.P. Sloan Foundation; the Bundesministerium 
fiir Bildung und Forschung, Germany; the Korean Sci¬ 
ence and Engineering Foundation and the Korean Re¬ 
search Foundation; the Particle Physics and Astronomy 
Research Council and the Royal Society, UK; the Russian 
Foundation for Basic Research; the Comision Interminis- 
terial de Ciencia y Tecnologia, Spain; in part by the Eu¬ 
ropean Community’s Human Potential Programme un¬ 
der contract HPRN-CT-2002-00292; and the Academy 
of Finland. 



FIG. 1: The measured inclusive jet differential cross section 
corrected to the parton level compared to the NLO pQCD 
prediction of the EKS calculation using CTEQ6.1M. 
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TABLE I: Results for the inclusive jet cross section corrected to the hadron level, d 2 <J hadron / dprdy , and to the parton level, 
d 2 a parton /dp T dy are shown for each bin together with the statistical (first) and systematic (second) errors. The correction 


factors, C h ^ p , applied to the hadron level 
additional 6% luminosity uncertainty. 

cross section to obtain the parton level cross section 

are also shown. There is an 
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